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Abstract 

We use factorized L operator to construct an integrable model with open bound- 
ary conditions. By taking trigonometic limit(r — ► y/—loo) and scaling limit(o; — > 0), 
we get a Hamiltonian of a classical integrable system. It shows that this integrable 
system is similar to those found by Calogero et al. 



PACS:75.10. 0530. 
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I Introdction 

In the last decades, a series of one-dimensional integrable many-body systems have been 
found [1-6]. One effective method to prove the integrability of the many-body systems is 
the Lax representation, which means that we can construct the complete set of integrals 
of the motion. The Lax representation for the elliptic Calogero-Moser model was found 
by Krichever [7]. And the Ruijsenaars-Macdonald's commuting difference operators can 
also prove the integrability of the many-body systems [5, 6, 8, 9]. 

It is known that there are a lot of two-dimensional exactly sovled models in statistical 
mechanics [10- 15]. The integrability is proved by the commuting transfer matrix. Recently, 
some relations between the one-dimensional many-body systems and the two-dimensional 
sovable models are found, see references [16-19] and the references therein. Hasegawa [16] 
found that the L operator [20-22] for the two-dimensional Z n Belavin model can be related 
to the Krichever's Lax matrix [7] . And the commuting difference operators given by the L 
operators is similar to the Ruijsenaars-Macdonald's difference operators. Hasegawa stud- 
ied the Z n Belavin model by imposing the periodic boundary conditions. We will study 
the Z n Belavin model by imposing the open boundary conditions, for the case of open 
boundary conditions see references [23,24] and the references therein. Using the factorized 
L operators and one solution of the reflection equation, we can constrcut a cummuting 
difference operators which should be equivalent to the Ruijsenaars-Macdonald's difference 
operators. By taking a special limit, we find a trigonometic integrable model which is 
similar to the trigonometric model found by Calogero et al [2-4]. Principally, we can 
obtain a classical integrable system if we can find a solution to the reflection equation. 
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The paper is organized as follows: In section 2, we will introduce the Z n symmetric 
Belavin model. The factorized L operator will be given in section 3. The commuting 
difference operators connected with the transfer matrix with open boundary conditions 
will be given in section 4. The special limit is taken in section 5, the integrable model is 
found in this section. Section 6 containes a summary and discussions. 

II The Z n symmetric Belavin R matrix 

The Z n symmetric Belavin R matrix [14,15] is given as 

Rl2 = ~ E W a { Z )I a ®I-\ (1) 

71 aezi 



with 



W a (z) = a A*±A , /q = g ^h a \h l3 = 5 i+ld , 



9a = u o i:j ,uj = e n e Z n ), 



a a (z) = 9 

t) = ^ e n ^^( m+a ^ T+27r ^~ I< - m+a ^ z+b ' ) (2) 
This R-matrix satisfy the Yang-Baxter equation(YBE) [10,11], 

#12 (zi - z 2 )Ri 3 (z 1 - z 3 )R 23 (z 2 - z 3 ) = R 23 (z 2 - z 3 )Ri 3 (zi - z^R^^ - z 2 ), (3) 

where R 12 (z) , R 13 (z) and R 23 (z) act in C n ® C n ® C n with R 12 (z) = R{z) <g> 1, R 23 (z) = 
1 <S> R(z), ect. One can find that the R-matrix also satisfy the following unitary and 
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2 n 

1 I "2 

2 ~ n 




cross-unitary properties. 



#12(21 - z 2 )R 2l (z 2 - zi) = p{z l - z 2 ) ■ id. 



(4) 



R 2 \(z 2 — zi— nw)R%(z 1 - z 2 ) = p(z 1 - z 2 ) ■ id. 



(5) 



where 



P(z) 



h(z + w)h(—z + w) 



(6) 



P{z) 



h 2 (w) 
h(z)h(—z — nw) 



(7) 



h 2 (w) 



w is defined by w = 77.77 an( i = ^(-z), ^ means the transposition in the i-th space. 
Assume an operator matrix L(z) satisfy the Yang-Baxter relation(YBR) 



with Li(zi) = L(zi) ® l,L 2 (z2) — 1 (g) L(Z 2 ). For the periodic boundary conditions, its 
transfer matrix is defined by = trL(z). We can prove that this tranfer matrices with 
different spectrum commute with each other [t(zi),t(z 2 )] = 0. For the open boundary 
conditions, Sklyanin [23] proposed a systematic approach to handle the problems which 
involves the reflection equation (RE) 

#12(21 - Z2)Ki(z 1 )R 2 i(zi + z 2 )K 2 (z 2 ) = K 2 (z 2 )R X2 (z x + z 2 )K 1 (z 1 )R 21 (z 2 - z t ), (9) 

where the reflecting K matrix is a solution of the RE. In order to construct the integrable 
models, we also need a dual reflection equation(DRE) which is associated with the cross- 
unitary relation of the R matrix. For Z n belavin model, the dual RE takes the form 



R 12 (zi - z 2 )L x {z^)L 2 {z 2 ) = L 2 (z 2 )L 1 (z 1 )R 12 (z 1 - z 2 ), 



(8) 



[24], 



#12(^1 - z 2 )K 1 (zi)R 21 (-z 1 - z 2 - nw)K 2 (z 2 ) 
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= K 2 (z 2 )R 12 (-Z! - z 2 - nw)K 1 (z 1 )R 21 (z 2 - z x ). (10) 

If we define the transfer matrix as t(z) = tr[K(z)L(z)K(z)L^ 1 (—z)], with the help of 
unitary and cross-unitary relations, one can prove that [t(zi),t(z 2 )] = 0, that means the 
model under consideration is integrable. 

Additionally, there is an isomorphism between K{z) and K{z) 

~ 777/; 

S:K{z)^K{z) = K{-z-—). (11) 

Given a solution K(z) of the RE, we can find a solution K(z) of the DRE. For Z n Belavin 
model, a solution K(z) to the RE is as follows [24], 

K(z)=K o (z)K o (0), (12) 

where 

K (z) = U 2a (z)u 2a ^I 2a , (13) 

°2o\,C) 

c is a arbitrary constant. 

Ill The factorized L matrix 

Jimbo et al [25] defined the intertwiner of Z n model as an n-element column vector 
4>a,a+fi(z) whose j-th element is 

4> ( <iLa(z) = 9 {j) (z-nwa,,nT), (15) 



6^\z — nwa^, nr) = 9 



1 _ L 

2 n 
1 

2 



[z — nwa^nr), (16) 



O// = CLfj, — j t Hu a v + Sp,a> £ Z n , S^s are some generic numbers. Using intertiwner, the 
face-vertex correspondence can be written as 



Rl2(zi - Z 2 )t/(p { a\-0,a-il(zi)(f>a- > f l ,a( Z ^ 



J2 W 



a — n — v a — jj, 
a — k a 



where W 
It is defined as 



a — ii — v a — /j, 
a — k a 



^-^^m^m, ( i7 ) 

z) is the face Boltzmann weight of A^li IRF model [25]. 



W 



a — 2fi a — jl 
a — fi a 



W 
W 



a — \x— v a — fi 
a — jl a 

a — fi — v a — v 
a — fi a 



(*) = 
(*) = 



h(z + w) 
h(w) ' 
h(z - <yw) 
h(-a^ u w) 
h(z)h((a flu + l)w) 



(18) 
(19) 



h(w)h(a IIU w) 

The other face Boltzmann weights are defined as zeroes. Where a^ v = — a v . At the 
same time, we can also find the n-element row vectors (f> and (f> which satisfy the following 
relations[20-22,25] 

®A>)*£U*) = <W> ( 2 °) 
^WZLM = <W (21) 



The above equation can also be written as 



Y J 4>a-H,a{z)4>a-a,a{ Z ) = J ' 
^2(j>a,a+t(z)(pa, a +t,(z) = I. 



(22) 
(23) 



Using those results, the face-vertex correspondence can be written in other forms 



(i — jl — k a — v 
a — ft a 



( Zl - z 2 )^!, k Jz 2 )^\^ a _,(z{), (24) 



a a + is 
a + ft a + fi + k 



{Z\ - Z2)(pa+iX,a+fi+k (*»), (25) 



if) 



a — n — v a — k 
a — fi a 



1 j 



[Zl ~ Z2)(pa2 k Jz 2 )^Jr l _ i> „ k ( y Z 1 ) 



(i') 



(26) 



u{Zl)4>ai+a{z2)Rl2{Zl - Z 2 ) 



5> 



a a + k 
a + fi a + fi + v 



{zi ~ ^)^i+U+/i+/>( 2 2)rf'a ) +«(^i)- 



(27) 



With the help of those face-vertex correspondence relations, we then can construct the 
L(z) matrix [20-22] which meet the YBR. Let 



f(a, fx, z){ = 4>%, a {z + £i)<^4 a (* + 6) 



(?) 



(28) 



£i and are arbitrary complex numbers. Then from the face-vertex correspondence 
relations Eq.(|TT|) and ( p6| ) , we have 



R(zi - z 2 )l- f(a - /}, is, z l Y- / 'f(a, is, z 2 ) 3 j , + £t < — ► v 



f(a - fa is, z 2 )] /(a, is, Zx)iR{z! - z 2 %?, + /i 



< > is, 



(29) 
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fi < — > v means the same form as the term before them, while the fj, and v exchange to 
each other. We introduce the difference operator T M 

r /t /(o) = /(a + /i)r /t , (30) 

and we have 

R(zi - z 2 )i-'T v f(a, u, z 1 )^V il f{a, u, z 2 % + \i < — ► v 
= r„/(a, v, z 2 y-T fl f(a, v, zi) % - Rfa - z 2 %'f + fi < — ► v. (31) 

We can find that the above equation is sufficient to prove the following YBR 

R( Zl - z 2 y; j j 'L(a,z 1 ) i ;L(a,z 2 y.; = L(a,z 2 ) j j 'L(a,z i y:R(z 1 - z 2 %'f. (32) 
Where L(a, z k ) = £ M r M /(a, n, z k ). Let 

g(a, z)\ = <P^(z + 6)0S +A (^ + 6) (33) 

There is 

£ r„/(a, !/, ^' ]T T_^(a, z)J = 5 ife (34) 

V fj, 

It shows that L _1 (a, z) = £ M r_ M #(a, z) . Then the transfer matrix with open boundary 



conditions is 



t(z) = trK(z)L(a,z)K(z)L~ 1 (a,-z) 

= trK(z) £ rV(a, u, z)K(z) £ r_^(a, /x, -z) 

V fj, 

= ^r^ivc^a,*), ( 35 ) 



with 

F${a, z) = 4>%_ &>a+(l (z + ^K(z)i^(-z + 6) (36) 

fJJfa *) = fiaUi-Z + mMU-O^ + 6) (37) 

G^(a,z) = F«(a,z)F^(a,z). (38) 

IV The calculation of tansfer matrix t(z) 

Substituting the if (z)(|l2D into Eq.(|36|) , we get 

= ^l^ a+fl (z) £ C/ 27 (z)a;^^[^ 27 ^ 271 ^o(0)]I^ +A (-^). (39) 

76^2 

Here we let £i = £2 = 0. We can prove that 

/-Vif o (O)0„, a+A Oz) = (-l) /3 - 1 e 2 -^^ +z -"^ + ^Va,a + A(-^+2^a, -ar-/?), (40) 

(see Apendix). The above Eq.(^) can be written as 

F$(a,z) = £ ^(z^^^W-^+i) 

x 0iV^(^)rf,U^ + 2 ^-27ir- 2 72 -l) (41) 

From Eq.(^) we know 4> a +ji-D,a+{i{z) can be obtained from the inverse of the matrix M 
whose elements are 





O^-^oa + ^-I)). (42) 



So we get 



detM' 

'a+A-^+A^C+A^ + Znwap - 2 7l r - (2 72 + 1)) = — ~r- (43) 



Substituting the v column elements of matrix M with the corresponding elements of 
column vector (p a ^ a+ p_(z + 2nwa fl — 2 7 ir — (2 7 2 + 1)) while keeping other matrix elements 
unchanged, we get the matrix M' . 

Suppose the elments of a matrix A defined as A$j = 9^'(nzj), one can prove that the 
determinant of the matrix A have the results [20-22] 

71 — 1 

detA = C{r)h(£ Zi - ——) ]J h( Zi - z k ). (44) 

i i<k 

Using this result, we get 

W^a.U^ + 2nW ^ ~ 2 ^ T ~ ( 2 ^2 + 1)) 

h(-z + wS + w(l-n) + ^- w{a v + <V - 1 + aj + ^(2 7 ir + 2 72 + 1)) 

h(-z + wS + w(l - n) + 2=i) 
x h(-w(a 3 + 5 M - 1 + <g + ^(2 7l r + 2 72 + 1)) 
h{-w(dj + 8^ -Ou- 5^)) 

(5 = J2i With the help of the formula 

a (z + -(ar + /?)) = e- 2 ^^ +z+ ^a a> p{z, r), (46) 



n 



finally, we get 



- W^27m 27r\/— T^2L 



o"2 7 i,2 7 2+i(-- 2; + w<5 + - n) + ^ - ^(^ + 8^-1 + a^)) 



x 



(7 (-z + w8 + w(l - n) + ^ 
a2 7l ,2y 2 +i(-w(dj + 8^ - 1 + a M )) 



x TT "^71^72 + 11 W \»J T » w .L-ru,^; 



Substituting the K(z) into Eq . (|37|) , considering the ismorphism relation between K(z) 
and -RT(^)(]ITD, we get 

F$(a, z) = ^(-z) £ U 2l {-z - ^)^^^[/ 7 ^ 27l ^o(O)]L0iVp,a +A (^) (48) 

7£^ 
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Using Eq fl40D , we have 



nW > 27 rV=T^p e 2 7 r^22L(^+a-n W (S,+V-l)+|) 



From Eq . (p3|) , we also know (j>a,,a+p,{— z) can be obtained from the inverse of the matrix 
M whose elements are 

= #(-z-moa A ). (50) 

So we get 

^al + a(-^ala-f>,a + u(-^ + 2 ™{a v + <W " 1) " 2 7i^ - 2 72 - 1) = (51) 



Substituting the \i column elements of matrix M with the corresponding elements of 
column vector 4> a +fi-p,a+fi(—z + 2nw(a u + 8^ — 1) — 2 7:L r — 2 72 — 1), holding other matrix 
elements unchanged, we get the matrix M' . Using this result of (fH|), we get 



C+At-^SA-^+A^ + 2nw & + ^ " 1) " 27i^ " 2 7 2 - 1) 
h(z + w5 + - w(q„ + - 1 + a„) + ^(2 7l r + 2 72 + 1)) 

/>(;? + u><5 + ^) 
T-r fe(-w(aj + - 1 + oj + ^(2 7l r + 2 72 + 1)) 
h{-w{aj - a M )) 

Similarly, with the help of Eq.(^Bj), we get 

= E U»y(-z - —) e - 2 ^—e 2 ^- 
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X 



o-2 7 i ,2 72 +i + w $ + r± 2 w(a ll + 5^ - 1 + a u )) 



x n 



a (z + wS+^) 

C r 2 7l ,2 72 +l(- W (Qj + ~ 1 + V)) 

cr (-w(aj - a,/)) 



(53) 



V Taking trigonometric and scaling limit 



When r — > v^~T°°) considering — 7^ m(m is integer number), we can find that 



1 , 2 



I m 

n 

1 1 272+I 



e 2^(^ + i+ m ), e 2.v^T22L 



Then, 



(54) 



(55) 



71 



sin7r(-z + u>5 + iu(l - n) + ^) 



x 1 1 - — t f = 

j+v smvr(a i + 5^ - a v - 5^) 



72 



(56) 



72 



72 



x ^ e' 27 "^ 1 " /™( r ) e 27r v / ^T(^ + i+m)(- 2 +n-i-n«,a M +272) 

^) 72 =0. 

We can see that from F^ 2 )(a,z), we can obtain the same result. So we only pay our 
attention to the case that ^ — m. The next step is to consider this case. 

We know there is an arbitrary parameter c in U2 7 (z). Let c = er + d with e < i, we 
have 



e 



then 



1 1 271 

2 ~ n 

1 , 272 

2 ~ n 



+ c, r) = e-W^^+s-K'+H 2 ^ 



1 + 271 +e 

2 ~ n ~ c 

1 I 272 

2 ~ n 



(z + c',r), (57) 



-27rV^Tez 



1 I 2 72 

2 "1" n 



(2 + c', r) 



1 , 272 

2 ~ n 



(C,r) 
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which have no relation to 72 and do not include any dynamics variables like a M , so we 
handle it as an irrelevant constant when taking trigonometic limit. The conclusion 
still hold when c = er + c'. Now, we can get(^- = m) 

sin7r(-z + w8 + w(l - n) + - w(a u + 8^-1 + a M ) + 



sin tt(—z + w8 + w(l - n) + ^] 



x -q smnj-wjaj + £ w - - 1 + cy) + gg) ^ 



sin n(-w(aj + 8^ — a v — 8^)) 



(2) x sin n(z + w8 + ^ - w(a^ + ^ - 1 + + 

nsin 7r(— w(a,- + £ U v — 1 + a u ) + 272+1 ) 
i V 3 . r-g s— (60) 
^ sm 7 r(-'u;(a i -a M )) 

For the n conservations do not depend on the spectrum parameter z, when z —>■ —y/—loo, 
we can get the integrable Hamiltonian of the system. So there are 



72 



x 



tt sin7r(-w(oj + <5 W - - 1 + oj + 232±1) 
sin 7r (-«;(% + - a„ - 8^)) 



72 

x tt sin7i(-w(fl J + ^ - 1 + cy) + g^tl) 
// M sin7r(-u;(a j - a„)) 

Expanding sin7i(— u>(a.,- + 8^ — 1 + a /t ) + 27 ^ +1 ) in F^)(a, z) and sin7r(— w(dj + 8^ — 1 + 

<^) + ~7r~) i n -^^( a i as the follows resepectively, 

' r e 7T V / ^(-«'(a J +5„-l+a M )+222±i) _ e _ 7rv ^T(_ u ,( 5j+ 5 Mj _i + a fl ) + H22±i)^ ^ 
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_} = ^V^l(-w(a j +S flu -l+a v )+^^-) _ e -7r v / ^T(-«>(a J +5 fl , y -l+a l/ )+222±i)^ 

sum over 72 and 7 2 and we can find that all the terms of Ff£)(a,z) and F[£)(a,z) have 
the forms 

Y^e 2 ^ lm " l \m = 1 ±1 ± ■ ■ ■ ± I = n, n - 2, n - 4, ■ ■ ■ , -n + 2, (65) 

and 

^e 2 ^ E!m7 2,m = 1 ±1 ± ■ ■ • ± 1 = w, w - 2, n - 4, • • • , -n + 2. (66) 

Only for case m = ±n, 0, the above summations are not equal to zeroes. When n = 
even,m can take 72, or 0, but it is rather complex to calculate the case m — 0. For 
simplicity, we only consider n = odd, so m can only take n. We have 

ffffo *) - e^^^^-i^+i) J] . . / jz * = t yT 5 (67) 

j>4 sinp2(-w;(a j + d Mj - o„ - 5^)) 

F<$(a, z) ~> e-W^EiH^+V-i+^H*) TT * (68) 

^sinn-^K- -°w) 

X TT 7 7^ r = T^T TT 7 7= "TT- ( 69 ) 

sin ir(-w(aj + d w - - a v - 6^)) ^ sin ir(-w(a k - a M )) 
For a^-afe = aj-afc+5,--4, and H = r_ /i r i/ G j[1I/ (a, z), in which r M /(a) = /(a+/})r M , 
we get 

1 



sin7r(-w;(a i - a„ + - <£„)) 



x ^[e-^^^)) J] — -7 7 ~ T~T ipd (70) 

^ sin 7r(-w;(a fc - a M + 5 fc - <y ) 



In quantum theory, one can set 



n d 



e A4/(x) = /(x + -L/i)ev^4. (71) 

V-l 
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ft d 

Compared with r M /(a) = f(a + fi)T^ let T M ~> e^ 9 *^, — V— l^tfy + ^) ^ ^ M and 
w ^ h, there is 



^sinhvr^ - x v ) 



x e ^ e 



sinh 7r(x fc - x^) ^' ^ 

When p\, it becomes a Hamiltionan of a classical integrable syetem. 

We also can begin directly from Eq.(|59|) and (|60f). Acorrding to the expanding method 
of getting the results fl67|) and (p8|), we have 

F^J(a.z) | e fv /r T(-2+ 2 2 1 + 1 - ,J ™f) - e -Tv /3 T(-^+V-+ 1 - n ^M)] 

x TT 7 7^ = ? j — tt, (73) 

f£ sm7r(-'«;(a i - a v + - d/xi/)) 



/if 



X 



n — it. (74) 



sin 7r(— u; (aj — a M )) 



while ignoring those sin functions which have no any dynamics variables. Let 

g>xv sin 7r(-to(aj + 5 Mj - a„ - 5^)) s i n ^(-^(a* - &„)) ' ^ 

we have 



(d ^) = |g 7r v^T( n + 1 - nu '( a /J+ a f +<W~1)) _|_ e -7r v / ^T(n,+l-niu(a fl +a^+(5 M y-l)) 

_ e 7r^ 3 T(-2,z-nu>(a, J +l-a„-<5 M „)) _ e -7i- v /3 T(-2z-n«j(a, I +l-a I ,-5 flJ ,)) |^ (76) 
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For different z, H{z) commute with each other, so we can obtain three commutable 
Hamiltonians H l from H 

W = Y J T~^uG^i = 1,2,3, (77) 

with 

Ql _ j e 7rv / ^(rt+l-mu(a M +a ! ,+5 M1 ,-l)) _|_ g-7r v / ^T(n+l-nty(a M +a, / +(5 M , / -l)) |^ (78) 
Q% _ e K\/^l(-2z-nw(a IJ ,+l-a v -8^)) g (79) 

In G^j,, if we have 5 M — > 8^ + p, where p is cummtable, we may still get two integrable 



Hamiltonians 



i/ 11 = ]Tr_,iv^, (81) 



H 12 = (82) 



with 

gfll _ g7r v /3 T(n+l-n«)(a M +a l/ +<5 Ml/ -l))^ (83) 
Ql2 _ e -ir\/^l(n+l-nw(a )1 +a l ,+S flv -l)) g (84:) 

H 11 and H 12 all can commute with H 2 and if 3 , but they may not commute with each 
other. 

An integrable system require its poisson bracket H — f±, fy, • ■ ■ , f n satisfy 
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When we take variables change, i.e. — > q\ = XQi + Pi> Pi remain unchange, in the new 
variable set {q^Pi}, the above equation fl85|) still hold, so the system is still integrable. 
Taking some kinds of limit after the variable subsitution, we may acheive a new integrable 
system, i.e. in Eq.([72|), let Xj — > x'^ = x x ji there is 



x e 



■It.' tt 1 



^ n -rrj-^-r,]- (86) 



sinh - (x' h — x 



When x ~ > 00 ) take the lowest order term, we obtain 

H" = Y e p »- p " TT TT . (87) 

This may be an integrable Hamiltonian. 

Additionally, we can let xq — > xo + po with others unchanged in H, and let sinh(7r(xo — 
Xi)) become le^ x °~ Xi ' to obtain a new Hamiltonian. We also can change two variables or 
more to get a new Hamiltonian. i.e. let 5q ^> Si 3> ■ ■ ■ 3> <5 n -i 3> 1 — >■ 00, Xi — > Xj + 5j, 
there are sinh7r(x; — Zj) — ► | e 7r ^ i ~ :!:j+<5l ~ 5j ^(i < j). Then 



x Y[ e~ n ( Xj ~ x ^ JJ e~^ Xv ~ x ^ JJ e~ n ( Xk ~ x " S) JJ g-^C^ -3 *). (88) 

j<V j>U k<[l fc>/i 

Like this, we can get amount of integrable Hamiltonians. 

VI Summary and discussions 

We study the Z n Belavin model with the open boundary conditions. By the factorized 
L operators, we constructed the commuting difference operators which is similar to the 
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Ruijsenaars-Macdonald's difference operators. We obtained a one-dimensional many-body 
integrable model from this commuting difference operators. By taking the special limit, 
this model is similar as the Calogero model. 

The one- dimensional integrable models presented in this paper is obtained from the 
transfer matrix operators with open boundary conditions. Generally for different solu- 
tions of the reflection equations, we can define different transfer matrices. Thus, to some 
extent, we give a systematic approach to obtain one-dimensional many-body integrable 
models. Besides the Z n Belavin model, there are a lot of two-dimensional exactly solvable 
model in statistical mechanics, it is interesting to find relations between those models 
with some one-dimensional many-body systems. 
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Appendix 



Using the formula 



[/^o(O)0 a>a+A (z)] w 

u^SijSj+^kdk+iflQ® ( z ~ nwa^, nr) 
[z — nwa^, nr) 



2 ' n 



(-l)u/V^^0 



1 

2 

1 i+a 

' 2 n 
_ 1 

2 



(— z + nwa^, nr) 

1 i+a 



2 n 
1 



-z + nwa^, nr). 



nr) 



i— a 
n 



(z,nr) (89) 



there is 



;-l)^- 1 ) e 27r ^ T ^(¥+^-" w ^+|) a; ^ 



2 n 
1 



(—2 + nwa^ — ar — (3, nr) 



(-1)<* 
(-1)<* 



-1) W=Tf^+ Z -n™a M +!)„i,,» 



-1) Jhr^TTafffli+^-rwua^+I) 



[fl$a,a+A( -JZ! + 2ra/;a / , - ar - (3, nr)] w (90) 



move the g in the right hand of the equation to the left side and we get the conclusion 
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